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ABSTRACT 

If ultra-high-energy cosmic rays (UHECRs) originate from extragalactic 
sources, understanding the propagation of charged particles through the mag- 
netized large scale structure (LSS) of the universe is crucial in the search for the 
astrophysical accelerators. Based on a novel model of the turbulence dynamo, we 
estimate the intergalactic magnetic fields (IGMFs) in cosmological simulations 
of the formation of the LSS. Under the premise that the sources of UHECRs are 
strongly associated with the LSS, we consider a model in which protons with 
E > 10^^ eV are injected by sources that represent active galactic nuclei located 
inside clusters of galaxies. With the model IGMFs, we then follow the trajectories 
of the protons, while taking into account the energy losses due to interactions 
with the cosmic background radiation. For observers located inside groups of 
galaxies like ours, about 70% and 35% of UHECR events above 60 EeV arrive 
within ~ 15° and ~ 5°, respectively, of the source position with time delays of 
less than ~ 10^ yr. This implies that the arrival direction of super-GZK pro- 
tons might exhibit a correlation with the distribution of cosmological sources on 
the sky. In this model, nearby sources (within 10 — 20 Mpc) should contribute 
significantly to the particle flux above ~ 10^° eV. 

Subject headings: cosmic rays - large scale structure of the universe - magnetic 
fields - methods: numerical 
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Introduction 



Over the past several decades, significant progress has been made on both theoretical and 
observational fronts in understanding the nature and origin of ultrahigh energy cosmic rays 
(UHECRs), those with -E > 1 EeV (= 10^* eV) (for recent reviews, see iNagano fc Watson 
2000l : iBerezinsky et a/.ll2006l ). Yet, the acceleration mechanism and the corresponding as- 
trophysical "accelerators" of these energetic particles are still largely unknown. Observa- 
tional data from several experiments such as the High Resolution Fly's Eye (Hi Res) indicate 



that the mass composition of UHECRs becomes lighter at higher energies (lAbbasi et al. 



20051 ) ■ However, composition analyses that include high energy interactions are often model- 
dependent and inconclusive (jWatsonll2006l ). According to a recent report from the Pierre 
Auger O bservatory, the m ass composition is likely mixed, possibly becoming heavier above 
30 EeV (jUnger et al. 2007 ). The ove rall distribution of UHECR arrival directions is consid- 
ered to be consistent with isotropy (IBurgett fc O'Malleyll2003l ). Exceptions to this general 
isotropy include the small-scale clusterings of doublets anc i triplets found in data collected by 



the A keno Giant Air Shower Array (AGASA) experiment (ITakeda et a/.lll999l : lUchihori et al 



2000). In additi on, a possible correlation o f AGASA events with BL Lacertae objects has 
been sugge sted (ITinvakoy fc Tkachevll200ll ). But these claims have been confirmed neithe r 
by HiRes (lAbbasi et all l2006l ) nor by Auger (IHarari et all 120071 : lArmengaud et all 120071 ) . 
However, the Auger Collaboration recently reported that the arrival directions of UHECRs 
above 60EeV in their data s how a correlation wit h the position of active galactic nuclei 
(AGNs) lying within 75 Mpc jAbraham et a/J 120071 ) . 



Since protons with E > 1 EeV cannot be confined within the Galactic plane, UHE- 
CRs likely originate from extragalactic sources. In particular, the overall isotropy of arrival 
directions suggest s that there may be a large number of sources dis tributed over cosmo- 
logical distances (INagano fc Watson I l2000l : iBurgett fc O'Malleyll2003l ). During their prop- 
agation through intergalactic space, such protons will lose energy by means of pion and 
pair production processes while i nteracting with the cosmic background radiation (IGreisen 



19661 : IZatsepin fc Kuzmin Ill966l ). The flux of ultra-high-energy (UHE) protons from cos- 



mological sources is thus expected to be strongly attenuated, resulting in a significant 
suppression in the obse rved spectrum above the GZK threshold energy, Eqzk ~ 40 EeV 
(IBerezinsky et a/.l 120061). Although the AGASA data show no indication of the GZK sup- 
pression (INagano &: Watson 1 12000| ). both the Yakutsk Extensive Air Shower array and 
HiR es have reported a s uppression of flux above Eqzk, contradicting the AGASA find- 
ing (jPravdin et all Il999l: IZech 1 120041) . Indeed the same suppression was seen in a recent 
Auger measurement (IFacal et a/.ll2007l ). which seems to have ended the controversy over the 
presence of the GZK c utoff. The so-called dip-c alibrated UHECR spectra from different ex- 
periments compiled by IBerezinsky et al\ (120061 ). appear to be in good agreement with each 
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other and to be consistent with GZK suppression. However, it has yet to be understood 
whether this suppression is actually due to the GZK cutoff or due to the maximum acceler- 
ation energy, i^max, of astrophysical accelerators. If UHECRs are protons, the GZK energy 
loss should opera te at acceleration sites as well, leading to an -Emax close to -Egzk (see, e.g., 
Kang et alhm^ . 



Most UHE protons observed above the GZK energy must come from within the so- 
called GZK sphere of radius -Rgzk ~ 100 Mpc, although the proton interactio n length at 
Eqzk is /40EeV ~ 1 Gpc, Corresponding to ^ ~ 0.2 (IBerezinsky fc Grigor'evalll988l ). However, 



finding cosmological sources inside the GZK sphere from the arrival directions of UHECRs is 



Siel et al. 


2003; 


Dolae et al. 


2004) 



(e.g. 

propagation of UHEC Rs in a structured an d magnetized universe, a dopting a numerical 
model for the IGMFs dSid et all [2003I . booi krmengaud et a/.ll2005h . In this model, the 
IGMFs are generated by means of the Biermann "battery" mechanism at shocks and then 



evolv ed passively in a cosmological hydrodynamic simulation (IKulsrud et a/.lll997l : iRyu et al. 



19981 ). The strength of the resulting fields is rescaled to match the simulated field strength in 
Coma-like clusters to the observed strength, which is on the order of microgauss. UHECRs 
with E > lOEeV are then injected at cosmological sources. These particles propagate 
through the magnetized large scale structure (LSS) of the universe, and arrive at a mock 
observer with deflection angle 6, the angle between the arrival direction and the source's 
location on the sky. Sigl et al. found that the deflection due to IGMFs is significant, with 
6 > 20° above 100 EeV. On the other hand, Dolag et al. adopted an IGMF model from a 
"constrained" cosmological simulation employing a magnetohydro dynamic (MHD) version of 
a smoothed particle hydrodynamics (SPH) code. They found the defiectio n angle of protons 
with 100 EeV to be less than 1°, contradicting the estimate by Sigl et al. (IDolag et a/.l 12004 
2005h . 



This controversy over the predicted deflection angles demonstrates the importance of 
modeling the IGMFs in identifying the astrophysical sources and studying the origin of 
UHECRs. In order to ree xamine thi s issue , here we adopt IGMFs, based on a novel models 
of the turbulence dynamo iRyu et al\ (120081 ) In this model, the strength of the magnetic fields 
is estimated from the local vorticity and turbulent kinetic energy in cosmological structure 
formation simulations. For the field direction, the passive fields from these simulations are 
used. 



The maximum energy of nuclei of charge Z that can be confined and accelerated by 
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astrophysical sources is given by 



-'-'max ~ 



-) 



L 
kpc 



lO^EeV, 



(1) 



where V, B, and L are the charact eristic flow speed, magnetic field strength, and linear 
size of the accelerator, respectively (iHillad Il984l ). There are a f ew viable candidates that 



100 EeV: jet s from AGNs fe.q.. iBiermann fc Strittmatter 



can p roduce the required -Emax 
1987), gamma-ray bursts {e.g., 
1997l) . In this study, we consider AGNs inside galaxy clusters as the sources of UHECRs, 



Waxmanlll995l ). and cosmological shocks (IKang et al. 



19961. 



Thus, the source position is in effect correlated with the LSS of the universe. Protons with 
E > lOEeV are injected at the sources and travel through the simulated magnetized space 
until they lose energy down to 10 EeV, meanwhile visiting mock observers placed inside 
groups of galaxies. In this study we focus mainly on the defiection angle, the time delay 
relative to rectilinear flight, and the energy spectrum of the UHE protons. 

In the next section, we describe our model IGMFs, cosmic-ray sources, and observers, 
and the simulations of the propagation of UHE protons in intergalactic space. The results 
are presented in §3. Finally, we conclude in §4. 



2. Models and Methods 

2.1. Intergalactic Magnetic Fields 

Our cosmological hydrodynamic simulations of a concordance ACDM universe are car- 
ried out with the following parameters: f^BM = 0.043, f^DM = 0.227, and Q\ = 0.73, 
h = ifo/(100 km/s/Mpc) = 0.7, and ag = 0.8. A cubic box of comoving size 100 h~^Mpc is 
simulated using 512^ grid zones for gas and gravity and 256^ particles for dark matter. In 
the simulation, magnetic fields are generated through the Biermann battery mechanism at 



structure formation s hocks and then evolved passively with the fiow motions (IKulsrud et al 



19971 : iRyu et a/.ll 19981 ). The simulation was repeated for six different realizations of the initial 
conditions to examine the effects of cosmic variance. Non gravitational effects, including ra- 
diative cooling, photo ionization and heating, and feedback from star formation are ignored. 
Those processes a. ffect the generation and evolution of magnetic fields mainly on small scales 



( iKang et a/.ll2007l ). which should not alter the large-scale fields primarily responsible for the 



deflection of UHECRs. 

In principle, if we were to perform full MHD simulations, we could follow the growth 
of the IGMFs through stretching, twisting, and folding of field lines, the process known 
as the turbulence dynamo. In practice, however, the computational resources currently 
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available do not allow high enough numerical resolution to reproduce the full development 
of MHD turbulence: since the numerical resistivity is larger than the physical resistivity by 
many orders of magnitude, the growth of the mag netic fields saturates before the dynamo 
action becomes fully operative (IKulsrud et a/.ll 19971 ). So. lSigl et al\ (l2003l . |2004| ) rescaled the 
strength of their passively evolved fields in the postprocessing analysis to match the observed 
field strength in clusters of galaxies. This rescaling, which hinges on the observed field 
strength in the intracluster medium, does not necessarily result in correct field strengths for 



filaments, sh eets, and voids. On the other hand, in the MHD SPH simulations of lDolag et al. 



(12004 . 12OO5I ). the initial field strength was adjusted to obtain a microgauss level in clusters 
of galaxies at the present epoch. They demonstrated that their simulated cluster fields are 
consistent with various observations, such as rotation measure profiles and the total radio 
powers of cluster halos. A lack of observations, however, prevents their simulated fields 
in low-density regions from being tested against the real magnetic fields in filaments and 
sheets. In the MHD SPH simulations, the fiow motions can be resolved reasonably well 
in high-density regions, where the smoothing length is sufficiently small, ensuring adequate 
growth of magnetic fields. However, turbulence may not be fully realized in the low-density 
regions, which have large smoothing lengths, so the field strength in filaments and sheets 
may be underestimated in such simulations. 



In this study, we take a new approach detailed in iRyu et al\ (120081 ). If one assumes 
that magnetic fields grow as a result of the turbulence dynamo, their energy density can be 
estimated from the eddy turnover number and turbulent energy density as follows : 



'^cddy 



^turb- 



(2) 



Here the eddy turnover time is defined as the reciprocal of the vorticity at driving scales, 
teddy = 1/i^driv (cu = V X w), and is the conversion factor from turbulent to magnetic energy, 
which is determined from high-resolution MHD simulations of turbulence. For our model 
IGMFs, the number of eddy turnovers is estimated as the age of the universe multiplied by 
the magnitude of the local vorticity, that is, tage^^- The local vorticity and turbulent energy 
density are calculated from the cosmological structure formation simulations. The energy 
density given by equation (2) fixes the strength of the IGMFs, so our model requires neither 
rescaling of the field strength nor adjustment of the initial fields. As in the work of Sigl et 
al., we assume that the topology of the IGMFs in the LSS can be represented statistically 
by the topology of the passive magnetic fields in the cosmological simulations. 

Figure 1 shows a two-dimensional slice of the magnetic field strength in our model at 
the present epoch. The IGMFs are structured and well correlated with the weblike cosmic 
distribution of matter. The strongest magnetic fields, with B > O.l/xG, are found inside and 
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around clusters, while the fields are weaker in filaments, sheets, and voids. Overall, there is 
a correlation between field strength and gas density, as can be seen in Figure 2{left). In the 
regions of galaxy clusters, with Pgas/ (Pgas) ^ 10^, we find (B) ~ 1 /iG. For typical filamentary 
regions, with Pgas/(Pgas) ~ 10, the field has (B) ~ 10~ ^ G. By coniparis on, the average field 
str ength in filanients is found to be (B) ~ 10"^ G bv Isid et al] booi ) and (B) ~ 10-^° G 
by iDolag et al\ (120051 ) . The right panel of Figure 2 shows the volume fraction, df/d\ogB 
{solid line), and its cumulative distribution, /(> B) {dotted line) and /(< B) {dot-dashed 
line). In our model, the volume f illing factor for B > 10^^ G is /(> lO^^G) ^ 0.01. By 



compa rison, /(> 10^*^ G) ^ 0.1 in lSid et al\ t004 \ and /(> 10^*^ G) ^ 10"^ in lDolag et al 
(120051 ). Hence, the field strength in filaments in our model lies in the middle of the values 
from the models used by these two groups, that is, lower than Sigl et a/.'s but higher than 
that of Dolag et a/.'s. 



2.2. AGNs as UHECR Sources 



In the framework of bottom-up acceleration models, AGNs are the most studied can- 
didate astrop hysical accelerators to p roduce cosmic-ray nuclei beyond the GZK energy (for 
a review, see iBerezinsky et a/.ll2006l ). As noted in §1, it has also been suggested that cos- 
mological shocks may accelerate nuclei of charge Z up to an -Emax of a few times 10^^ Z eV, 



(Kane et al. 


1997; 


Inoue et al. 


2007) 



cal shocks can provide the acceleration sites for super-GZK cosmic rays (llnoue et a/.l 120071 ). 
In this paper, we consider only AGN-like objects as sources of UHECRs, since we focus on 
the propagation of UHE protons. 

We identify X-ray clusters characterized by the X-ray emission weighted gas temperature 
kT > 1 keV in the cosmological simulations as source locations. One AGN is placed at the 
center of each "host" cluster. With this selection criterion, the source locations are chosen, 
in effect, at high-density regions with the strongest field strength. Table 1 shows the number 
of such clusters found in the six cosmological simulations (ACDM1-ACDM6) with different 
initial conditions. The three-dimensional distribution of the 18 sources in ACDMl is shown 
in Figure 3. Given the simulation volume of (100 h~^Mpc)^, the mean separation of sources 
is Is ~ 40 h~^Mpc and the source number density is = 2 — 3 x 10~^/i^Mpc~^, which is 
consistent with the require d UHECR source density inferred frorn the small-scale clusteriri g 
found in the AGASA data JVoshiguchi et a/.ll2003l : ISid et all^W^ . iBlasi fc De Marcoll2004[ ). 
The field strength at the source locations mostly lies in the range 0.1/xG ^ Bg ^ 2/xG, with 
a peak at ~ l/xG; its distribution is shown in Figure 4 {left). 
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2.3. Groups of Galaxies as Mock Observers 



The key physical condition for an "observer" that is most relevant to this study is 
the strength and direction of the magnetic fields, since we are interested primarily in the 
defiection angles and time delays of UHECRs. Little is known about the magnetic fields in 
the intergalactic space within the Local Group. So, we select groups of galaxies identified 
in the simulation data that h ave similar halo gas teniperatu res to the Local Group, that 



are located in magnetic environment similar to that of the Local Group. There are about 
1000 — 1400 identified groups with gas temperatures in this range inside the simulation 
volume (see Table 1). As can be seen in Figure 3, these groups are not distributed uniformly 
but are located mostly along filaments, following the matter distribution of the LSS. A mock 
observer, modeled as a sphere of radius i?obs = 0.5 h~^Mpc, is placed at each group. The 
value of i?obs is chosen so that the observer's sphere is well contained within the associated 
filament, since the typical thickness of the magnetized region around a filament is about 
2 — 3 h^^Mpc. If we were to use an i?obs smaller than the value adopted, smaller cross 
sections would lead to smaller detection rates of cosmic rays in our numerical experiment 
described below. The right panel of Figure 4 shows the field strength within the observer 
spheres for the six cosmological simulations: 10~^/iG < -Bobs ^ 0.1/iG with a peak at ~ 10~^ 
G. This illustrates the distribution of magnetic field strength in filaments. 



The characteristics of the particle spectrum accelerated at UHECR sources are largely 
unknown. We thus make the simple assumption that the accelerated protons have a power- 
law energy spectrum at their source: iVinj(-Einj) oc E^^ for 10 EeV < E^^^ < 10^ EeV. In 
practice, 7 = is used to generate a fiat injection spectrum at the sources, and later a 
weighting factor proportional to E^J is applied to the statistics (except for the distribution 
in Fig. 5; see §3.1for details). For each ACDM simulation, a total of 3 x 10^ particles are 
randomly distributed over the sources and then launched in random directions from random 
positions inside a sphere of radius 0.5 h~^Mpc. 

We follow the trajectories of the UHE protons by numerically integrating the equations 
of motion in our model IGMFs, 



is, 0.05keV 




assuming that these groups 



2.4. Propagation of UHE protons in the Intergalactic Space 




(3) 



The energy losses due to photo pion production and pair production are treated with the 
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continuous loss approximation (IBerezinsky et a/.ll2006l ). but the adiabatic losses due to cos- 
mic expansion are ignored, because the largest source-to-observer distance is -Dmax ~ 1 Gpc 
in our experiment, corresponding to only ^max ~ 0.2. In practice, the distances are not 
known in advance, since we are integrating the trajectory from sources to observers. 

The simulation box of (100 /i~^Mpc)^ at 2; = is used to define sources (host clusters), 
mock observers (groups), and the IGMF data. Additional virtual boxes with the same distri- 
bution of mock observers and IGMF data are periodically stacked indefinitely. Particles are 
injected from the sources only in the original box. Then they travel through the magnetized 
space consisting of the original box and the rephcated periodic boxes. Once a particle visits 
an observer sphere, the arrival direction, time delay and energy of the particle are registered 
as a "recorded event." We let the particle continue its journey, visiting several observers 
during its full flight, until its energy falls to 10 EeV. With 3 x 10^ protons injected for each 
simulation box, about 2.6 x 10^ events are recorded in total for all six ACDM simulations. 

In our propagation experiment, the source-to-observer distance D can be arbitrarily 
small because the specific way that we set up source and observer locations, and the statistics 
of recorded events depend on t he minimum value o f D. In reality, the closest AGNs to us 
are Centauru s A, at 3.42 Mp c (jPerrarese et a/.l 120071 ). in the southern hemisphere and M87, 
at 16.7 Mpc (IMei et a/.ll2007l ). in the northern hemisphere. So we mostly present the results 
for -Dmin = 3 Mpc or ii'min = 10 Mpc. The recorded events with D < D^i^ are excluded from 
the analysis. For -Dmin > 10 Mpc, however, the results become less sensitive to the value of 
D 

min • 



3. Results 

3.1. Deflection Angle and Time Delay of UHECRs 

With a gyroradius 

UHE protons will suffer significant deflection during their propagation when they pass regions 
with B > 10^® G, that is, clusters and filaments. Here the filamentary regions are more 
significant players than clusters, since the volume filling factor of filaments is much larger 
than that of clusters (see Fig. 2). As a result of this deflection, the actual path traveled 
by the UHECRs in the presence of the IGMFs can be much longer than a rectilinear path, 
causing a significant time delay. We therefore measure the deflection angle as the angle 
between the arrival direction of the cosmic rays and the source position on the sky, and the 
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time delay as the difference between the arrival time and the rectilinear travel time. 

UHE protons can be also deflected inside the host clusters of sources, before they escape 
to intergalactic space. Typical clusters have a magnetized core envelope structure with 
Score ~ 1a*G, -Rcore ^ 0.5 h~^Mpc and Senv ~ 0.01 — 0.1//G, extending out to it!env ~ 
3 h'^Mpc (see Fig. 1). So, the protons with E < Eqzk injected by AGNs are scattered 
by tiirbulent magnetic fields inside the host clusters and confined within the magnetized 
structure for a while. The scatterings by the turbulent fields local to the sources alone can 
cause a deflection angle 



for protons with E < Eqik- 

Figure 5 shows the distributions of the deflection angle 9 and the time delay td as 
functions of the distance D for the ACDMl simulation. The events are divided into three 
channels in observed energy, as follows: 10 < E'obs < 30 {top), 30 < E'obs < 60 {middle), 
and £^obs > 60 {bottom), where the particle energy is given in units of EeV. The data 
points are color-coded by the injection energy in the same three channels, that is, red for 
10 < E'inj < 30, blue for 30 < E'inj < 60, and green for E'inj > 60. For these plots, 
the calculation was performed with an injection spectrum with 7 = 2.7 instead of the fiat 
spectrum, since the recorded data points in this type of representation cannot be weighted 
with a factor proportional to E^J . 

Sub-GZK protons, with Eobs < 60 EeV, come from sources as distant as D ~ 1 Gpc. 
For these particles, the distribution of 6 shows a pattern roughly in accord with the diffusive 
transport limit, but it also indicates a bimodality divided at D ~ 15 Mpc (Fig. 5, top and 
middle left). The events with D < 15 Mpc are likely to be cases in which both the source 
and observer belong to the same filament. These particles are more likely to travel through 
strongly magnetized filaments rather than void regions, resulting in large deflection angles 
in addition to a large 6'source given by equation (5). On the other hand, for the events with 
D > 15 Mpc the particles come from distant sources associated with different filaments. 
Some of these may fly through voids and arrive with small 9, while most are deflected 
significantly by the IGMFs. On the other hand, on average the time delay tends to increase 
with distance D, as expected. 

Super-GZK protons, with E^hs > 60 EeV come mostly from within D ~ 100 Mpc. Most 
of them from D > 15 Mpc arrive with 9 < 10° and t^ < 10^ yr (Fig. 5, bottom). Since the 
volume filling factor for B greater than a few times 10~^G (corresponding to ~ 2 — 3 Mpc) 
is small, those particles could travel almost rectilinearly through void regions, avoiding the 
strongly magnetized regions of clusters and groups. 




(5) 
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in Figure 6, the distribution of recorded events is shown in the [Eohs, d) and (-Eobs, logt^) 
planes. The events recorded in all six ACDM simulations are included and 7 = 2.7 and 
-Dmin = 3 Mpc are used. On average the deflection angle decreases with energy, but a 
clear transition from the diffusive transport regime to the rectilinear propagation regime 
is apparent around E'obs ~ -Egzk- Sub-GZK particles, with long interaction lengths and 
small gyroradii, are strongly scattered, while super-GZK particles, with short interaction 
lengths and large gyroradii, are much less affected. As expected, the time delay decreases 
with increasing energy on average but has a rather wide spread at any given energy. For 
comparison, the rectilinear flight time for the mean separation of sources. Is = 40 Mpc, is 
tree = Is/c ^ 10^ yr. 

Figure 7 shows the fractions of recorded events in all six ACDM simulations as functions 
of 9 and t^, df /dO and df /dlogtd, and their cumulative distributions. The events are divided 
into three energy channels as in Figure 5, and each curve is normalized by the total number of 
events in the corresponding channel. In order to demonstrate the dispersion of the statistics 
due to cosmic variance, we also plot the error bars in the cumulative distributions, which 
are calculated as the standard deviations of the values of / for the six simulations. 

In the lowest energy channel (10 EeV < i^obs < 30 EeV, red lines), the deflection angle 
is quite large, with about 70% of the events arriving with 9 > 30°, that is, /(> 30°) ~ 0.7. 
Moreover, with /(> 10® yr) ^ 0.7 the time delay is much longer than the typical lifetimes 
of AGNs {tagn = 0.01 — 0.1 Gyr). In the highest energy channel (Eobs > 60 EeV, black 
lines), on the other hand, about 70 % of the recorded events arrive with a deflection angle 
smaller than ~ 15° and a time delay less than ~ 10^ yr. About 35 % arrive with an angle 
smaller than ~ 5°. This implies that the arrival direction of super-GZK cosmic rays may 
show a positional correlation with the source AGNs and also with the LSS, and the source 
AGNs are very likely to still be active for such events. We note, however, that these results 
are not restricted to the speciflc AGN model. They can be applied for any UHECR sources 
that have a spatial distribution and magnetic fleld environment similar to those of luminous 
X-ray clusters. 

We note t hat the present work produces results difl Ferent from what previous studies 



have predicted (lArmengaud et a/.ll2005l : iDolag et a/.ll2005l ). Speciflcally, the deflection angle 



is smaller than that found by Sigl et a/.but larger than that of Dolag et ai. This should be 
attributable to the difference in the models for the IGMFs, as discussed in §2.1. We also 
note that the effects o f Galactic magnetic flelds are not included in our analysis. Recently, 



Takami fc Satol (120071 ). for instance, considered several different models for the Galactic 
magnetic flelds and predicted that the deflection angle of 10^^'® eV protons should be greater 
than 8° toward the Galactic center, while being mostly 3 — 5° outside a circular region of 
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30° radius around the Galactic center. 



3.2. Predicted Energy Spectrum 



Here we present the energy flux, J(i?obs), of the recorded cosmic-ray events in our 
propagation experiment. By applying a weighting factor proportional to E^J to each 
recorded particle, the energy spectra for different values of 7 can be constructed. The 
predicted spectra are calculated for injection spectra with 7-values of 2.0, 2.4, and 2.7 and 
for -Dmin = 3 and -Dmin = 10 Mpc. Again, all data from the six ACDM simulations are 
combined. Figure 8 shows the re sulting sp ectra, J(E), alorig wit h the data observed a t 



AGASA (INagano fc Watson 



2000), HiRes-I (IBerezinskv et a/.ll2006f ). HiRes-II fIZech Il2004h 



and Auger (jParizot et a/.l 120071 ). Since the amplitude of the injection spectrum is not spec- 
ified, the amplitude of the predicted J{E) is arbitrary in our model. Therefore each curve 
was adjusted by eye to fit the HiRes data below Eqzk- 

The presence of GZK suppression above 60 EeV is obvious in all the predicted spectra 
and the observed data except for the AGASA data. The predicted spectra for 7 = 2.4 — 2.7 
are all consistent with the observed data, again with the exception of the AGASA data. 

As shown in the bottom panel of Figure 8, our work predicts that above ~ 100 EeV 
the flux is much higher with -Dmin = 3 Mpc than with -Dmin = 10 Mpc, indicating that the 
contribution from nearby sources is important. Thus, if the injection spectrum has a power- 
law distribution extending well beyond the GZK energy as we assume here, the implication is 
that the Auger experiment, which has Centaurus A in its field of view, may see a higher flux 
of super-GZK cosmic rays, comp ared with experirnents in the Northern Hemisphere such as 
HiRes and the Telescope Array (iFukushima et a/.l 120071 ). However, it is quite possible that 
the injection spectrum is limited to a maximum energy -Emax set by the age and size of the 
astrophysical accelerators or by the GZK energy loss at the acceleration sites. Moreover, the 
value of -Emax may vary with the properties of the accelerator, rather than maintaining a 
constant value of 10^ EeV. In order to settle this issue, much better statistics for the energy 
spectrum and the arrival directions above 100 EeV are needed. 



4. Summary and Discussion 

In the search for the astrophysical sources of UHECRs, it is important to understand 
how the propagation of these charged particles is affected by intergalactic magnetic fields 
in the large-scale structure of the universe. On the other hand, the information imprinted 
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on the distribution of the UHECR arrival directions may help us to understand the nature 
of the IGMFs and their roles in the formation and evolution of the LSS and constituent 
galaxies. Considering the limitations of current observational techniques in measuring the 
IGMFs in very low density regions such as filaments and voids, it is crucial to construct a 
physically motivated model to estimate the IGMFs in the LSS. 



In this study, we adopted a new model based on the turbulence dynamo (IRyu et al. 



20081 ) to predict the strength of the IGMFs. The magnetic field energy is estimated from 
the local vorticity and turbulent kinetic energy of flow motions in cosmological simulations 
of the LSS formation in a concordance ACDM universe. For the direction of the IGMFs, the 
topology of the passive magnetic fields followed in the cosmological simulations is used. This 
approach provides an IGMF model that is independent of the initial seed fields and does not 
require any renormalization to yield the observed field strength in the intracluster medium. 
We predict highly structured IGMFs with characteristic field strengths on the order of 10~^ 
G in clusters of galaxies and 10~® G in filaments. The fields should be much weaker in sheets 
and voids. 

Protons with 10 EeV < Ei^^ < 10^ EeV are injected at the locations of luminous X-ray 
clusters with kT > 1 keV. These sources may represent a population of AGNs residing in- 
side host clusters. This X-ray temperature criterion naturally places the sources at strongly 
magnetized regions with Bg ~ O.l/xG with a comoving density of 2 — 3 x 10~^( h''^Mpc)~^ . 
Then the propagation of the UHE protons is followed through the structured IGMFs, in- 
cluding the energy losses due to interactions with the cosmic background radiation The UHE 
protons are recorded at the positions of mock observers located in groups of galaxies, with 
halo temperatures in the range O.OSkeV < kT < O.SkeV. 

Below the GZK energy, the UHE protons come from sources as distant as ~ 1 Gpc. 
They are significantly scattered by the IGMFs, resulting in a wide range of deflection angles, 
up to 180°, and have time delays ranging from 10^ to 10^'^ yr. On the other hand, the 
protons above 60 EeV come mostly from sources within ~ 100 Mpc. About 70% of them 
avoid strong deflection and arrive at the observers within ~ 15° of their source position on 
the sky with a time delay of less than ~ 10^ yr. About 35% arrive within ~ 5°. This implies 
that there may exist a correlation between the arrival direction of super-GZK cosmic rays 
and the sky position of the corresponding AGNs. We thus conclude that in the present 
scenario, UHECR astronomy may be possible at > 60 E eV. Our prediction se ems to be 



consistent with a recent report by the Auger Collaboration (lAbraham et a/.ll2007l ) in which 
the arrival directions of cosmic rays above 60 EeV in their data were found to be correlated 
with the sky position of AGNs within 75 Mpc. 

For any cosmological sources, we expect to see GZK suppression in the energy spectrum 
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of UHECRs if the injection spectrum has a power-law distribution and extends well beyond 
the GZK energy. In this case, nearby sources, within 10 — 20 Mpc, are expected to make a 
significant contribution to the flux above ~ 100 EeV. 



Finally, as recently reported by Auger (jUnger et a/.l 120071 ). some UHECRs might be 



heavy nuclei. In the cosmological-shock model, for example, protons can be accelerated up 



(Kane: et al. 


1997; 


Inoue et al. 


2007) 



of heavy nuclei from cosmological sources in our model IGMFs, taking into account photo 
disintegration, photo pair production, and photo pion production processes. We expect that 
the propagation of UHE iron nuclei [Z = 26), at least, will be i n the diffusive transpor t 



regime as a results of their much smaller gyroradius (r^ (x E/Z) (lArmengaud et a/.ll2005l ). 



For intermediate-mass nuclei such as He, C, N, and O, detailed propagation simulations 
including secondary particles produced by photo disintegrations are necessary in order to 
determine whether astronomy with UHE nuclei is possible or not. 
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Table 1. Numbers of Sources and Observers 



Simulation 


Sources 


Observers 


ACDMl 


18 


1000 


ACDM2 


31 


1344 


ACDM3 


20 


1379 


ACDM4 


29 


1336 


ACDM5 


24 


1343 


ACDM6 


28 


1365 
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Fig. 2. — Le/t: volume fraction in the gas density vs. field strength plane with our model 
IGMFs at redshift z = 0. Right: volume fraction, df/dlogB {solid line), and its cumulative 
distributions, /(> B) {dotted line) and /(< B) {dot-dashed line) as a function of the IGMF 
strength. 
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Fig. 3. — Distribution of UHECR sources {red circles) and mock observers {blue circles) 
in simulation ACDMl. Sources are modeled as AGNs inside X-ray clusters, while mock 
observers are placed inside groups of galaxies similar to the Local Group. 
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Fig. 4. — Left: Distribution of magnetic field strength at source locations in the six simula- 
tions with different initial conditions. Different line styles are used for simulations ACDMl 
- ACDM6. Right: Distribution of magnetic field strength within the observer spheres. The 
same line styles as in the left panel are used. 




Fig. 5. — Deflection angle (6) and time delay (t^) as functions of the source-to-observer 
distance (D). The events recorded with observed energies 10 EeV < Eohs < 30 EeV are 
shown in the top panels, 30 EeV < Eohs < 60 EeV in the middle panels, and i^obs > 60 EeV 
in the bottom panels. The data points are color-coded by the injection energy as follows: 
red, 10 EeV < E^^j < 30 EeV; blue, 30 EeV < E^^^ < 60 EeV; green, Ei^j > 60 EeV. An 
injection spectrum N{Eia}) oc E^^''^ is assumed. (See text for details.) 
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Fig. 6. — Distribution of observed UHECR events in the planes of observed energy vs. de- 
flection angle [left) and time delay (right). The events recorded in all six ACDM simulations 
are included, and 7 = 2.7 and -Dmin = 3 Mpc were used. 
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Fig. 7. — Fractions of observed UHECR events as a function of deflection angle {top left) 
and time delay {top right), and the respective cumulative distributions {bottom). The events 
recorded in all six ACDM simulations are included, and 7 = 2.7 and D^[^ = 3 Mpc were 
used. The distributions for different energy channels are shown: red, 10 EeV < Eohs < 30 
EeV; blue, 30 EeV < i^obs < 60 EeV; black, Eohs > 60 EeV. The error bars shown for the 
cumulative distributions are the standard deviations of / for the six simulations. 
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Fig. 8. — Energy spectra of UHE protons predicted by our model. The injection spectrum 
at the sources is proportional to E^^ for 10 < i^inj < 10^ EeV. Top: The blue, red, and black 
lines are for 7-values of 2.0, 2.4, and 2.7, respectively, with a minimum source-to-observer 
distance -Dmin = 3 Mpc. Bottom: The red lines are for an injection spectrum with 7 = 2.4, 
while the blue lines are for 7 = 2.7. The solid hues are for -Dmi n = 3 Mpc, and the dotte d 
lines are for -Dmin = 10 Mpc. The data observed at AGASA (INagano fc Watson 1 120001 ). 
HiRes-I JBerezinskv et allboOfil ). HiRes-II Jzech II2OO4I ). and Auger jParizot et oi.ll2007h are 
marked with asterisks, open circles, filled circles, and triangles, respectively. The predicted 
spectra were arbitrarily scaled by eye to fit the HiRes data below Egzk- 



